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1. Introduction 

Observables related to meson oscillations represent among the best measured quantities in 
particle physics. This is true particularly for the K — K and the B ( i — B,i cases, where experiments 
at K- and 5-factories have pushed precisions to the permil level and new quantities in the D and 
B s sectors promise to add up thanks to now-running and forthcoming experiments. Experimen¬ 
tal commitment is related to the special interest meson oscillations have on the theoretical side: 
they constitute eminent sources of so-called flavor changing neutral currents (FCNC’s) possibly 
accompanied by a sizeable violation of the CP symmetry. FCNC’s and CP violation are in turn 
the effects that most naturally (and copiously) occur in basically all extensions of the SM, and in 
particular in SUSY, that, notwithstanding the success of the SM, represents the “standard way” 
beyond it [[!]]. The most popular low-energy realization of SUSY is provided by the MSSM, that 
however features ‘by default’ a huge space of parameters. Placing high precision constraints on the 
latter thus becomes of capital importance, considering that it can provide key handles on possibly 
ungrasped symmetries. To this aim, meson oscillations constitute privileged observables, recalling 
their experimental and theoretical ‘virtues’, as mentioned above. 

The amplitudes for meson oscillations can be accessed in the Effective Hamiltonian (EH) 
approach, through the so-called A F = 2 EH’s, reading as 

^r =2 =L<w&, (U) 

i 

where the Qj are local operators of increasing mass dimension and the C, the corresponding ‘cou¬ 
plings’ (Wilson coefficients). In the MSSM and restricting to the A B = 2 case (B,j — B c f) a suitable 
basis for the Qj up to dimension 6 is given by []|] (Pr.l = (1 ± 7s)/2) 

Q\ = (WbtfVd) (WbYtiLVd), Qi = (WiAVd) (Wl p LVd ), & = OfifrVd) iWiPLWd), 
<24 = (Wb P L¥d) C Wl P RV J J ) , Q5 = C % P LVd ) iWb P RVd) , ( 1 . 2 ) 

plus Q 12,3 obtained from <2 i, 2,3 with the substitution L —* R. In eq. ( |P| ) indices i, j denote color. 

The Cj in eq. ( |1.1[ ) are calculated by requiring that the amplitude sf for meson oscillation, 
evaluated in the full (i.e. in the MSSM) and in the effective theory ( |1.1[ ) be the same below a high- 
energy threshold Msusy- Since the C,- encode the short-distance ‘part’ of the process, the equality 
4^fuii = ■ (/ Y<i it can be imposed by choosing external quark states with arbitrary kinematics. 

The C,(/i ~ Msusy) must then be evolved through RG flow to a scale pi ~ few GeV. The 
anomalous dimension for the operator basis ( |l.2[ ) needed in the RG equations is known to NLO 
accuracy [jj]. The amplitude for B c i — B ( j mixing is finally obtained by multiplying the Cffi) by the 
matrix elements ( Qi){p ) between the physical states of interest, evaluated on the lattice [4j]. 

The leading order (LO) C, for AF = 2 Hamiltonians in the strong sector of the MSSM were first 
calculated in refs. [Q]. They are simply obtained by computing four gluino-squark ‘box’ diagrams 
(analogous to the VT-quark boxes in the SM), and rewriting their amplitude in terms of tree-level 
matrix elements of the operators (|1 ,2|). 

It is clear that knowledge of the next-to-leading order (NLO) corrections to the C, is the missing 
ingredient for a full NLO analysis of AF = 2 processes. Such corrections are the subject of the 
present contribution. In particular, section [2] briefly presents the highlights of the calculation, that 
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is dwelt on in a forthcoming publication. Section [| addresses the phenomenological analysis, 
presenting preliminary results for the B d — B d system. A thorough analysis is again postponed to 
an upcoming paper. 

2. NLO corrections to A F — 2 Hamiltonians 

NLO corrections to the C, entering eq. 0 are calculated by enforcing the equality sffa\\ = 

at the scale Msusy up to two loops in the full theory and correspondingly one loop in the ef¬ 
fective one. The two-loop diagrams mainly consist of gluon or squark corrections to the LO boxes: 
an inventory taking into account multiplicities and ‘equivalences’ results in around 70 diagrams. 
Hence, to be on the safe side as for the generation and the initial treatment of the full amplitude, 
the Mathematica package FeynArts [|j] was an essential tool. 

To two-loop order it becomes important to choose a regularization scheme for ultra-violet 
(UV) divergences as well as a kinematical configuration for the external quark legs. Concerning 
the first issue, we performed the whole calculation in both the NDR and DRED (MS) schemes. 
As for external legs, we chose massless and zero-momentum quarks, since this leads to enormous 
simplifications in the algebra and in the treatment of the two-loop integrals, that in this way are 
just “vacuum" ones. This choice brings about infra-red (IR) divergences that however must cancel 
in the matching. In the intermediate steps IR divergences were coped with either by dimensional 
regularization (like UV ones) or by endowing gluon propagators with a mass term. 

We then compared the calculations in all the resulting UV-IR regularization schemes and ob¬ 
tained perfect consistency, taking into account the findings of [^] for the treatment of evanescents 
and of [j7p for the fact that NDR breaks SUSY. Other important checks of the calculated C, were the 
already mentioned disappearance of IR divergences and the fact that the /./-dependence fulfills the 
relevant RG equations. Explicit formulae will be reported elsewhere. 

We note that the final results are in the so-called mass insertion approximation (MIA) for the 
entries of the squark mass matrices in the ‘super-CKM’ basis: within the MIA diagonal entries 
are approximated with a common mass term A/ 2 appropriately chosen (1 st approx.) and the off- 
diagonal ones A xy (Y, Y being the four possible chiralities) are considered small with respect to M~ 
and expanded as mass insertions 8xy = A xy/M; <C 1 (2 nd approx.). The MIA enormously simpli¬ 
fies the phenomenological analysis and in our NLO case also decides for the analytical feasibility 
of the two-loop integrals. 

3. Phenomenological analysis of the B d — B d system 

The calculated NLO corrections to the C, at the matching scale, evolved to the scale /./ with the 
NLO ADM as reported in [||], permit for the first time a full NLO control on A F = 2 Hamiltonians 
(0). This allows a correspondingly accurate analysis of meson-antimeson oscillations, updating 
the results of [$]. We restrict here to the B d — B d sector and use the findings of ref. 0 for the 
lattice determination of the matrix elements for the effective operators ( |L2| ) in order to access the 
relevant amplitude stf Bd = (B d \(J4?Jf =2 ) SM+slJSY \B d ). 

Our analysis proceeds as follows: one extracts with normal distributions around their central 
experimental values all the (SM) parameters needed in sf Bi (CKM entries are taken from tree- 
level processes only, to leave SUSY unconstrained). The gluino mass is set to Mg = 500 GeV 
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and x = M^/M~ = 1, while the mass insertions are extracted with flat distributions according to 
the intervals: Abs(5xy) £ [0, 1] and Arg(Sxy) £ (-71, 7i\. Therewith one has MonteCarlo deter¬ 
minations for tl . that are compared with experiments via the observables AM r j = 2Abs(.c//y,,) 
and 2/3eir = Arg (g^B d ), the latter being accessed by studying the decays B t \ —> Jj\\iK s '. One can 
then place constraints on the mass insertions by ‘switching on’ one 8xy at a time according to the 
choices: 'LL only’, ‘RR only’, 'LL = RR\ ‘LR only’, ‘RL only’, ‘LR = RL'. Barring accidental 
cancellations among amplitudes proportional to different 8xy parameters, the results can be inter¬ 
preted as maximum ranges allowed by the present experimental information to the various mass 
insertions. Constraints on 8ll are displayed in fig. [j], left panel, whereas the right panel nicely 
shows the reduction in the scale dependence of .e/[j r/ when only LL and RR insertions are kept. 

The ‘hierarchical’ pattern of the constraints, e.g. max(t)/jj A> max(5«L), represents an a pos¬ 
teriori justification for the MIA, wherein possibly large interference effects among squark mass 
matrix parameters are not taken into account. Such pattern also motivates our separate considera¬ 
tion of LL, RR and respectively LR. RL insertions in the analysis of the scale dependence. 


-P-15 



Re((5 d , 3 ) LL 



Figure 1: Left panel: allowed ranges for Re(d/2) vs. Im(5^). The sides of the ‘boxes’ are proportional 
to the number of ‘events’ as weighted by the MonteCarlo. Constraints for different observables are super¬ 
imposed with different colors: AM c / only (magenta), sin(2j3) only (green), sin(2/3)&cos(2j3) (cyan), all the 
constraints (blue). Right panel: scale dependence of the LL and/or RR terms in the amplitude in the cases of 
Wilson coefficients to LO or to NLO accuracy (similar results hold in the LR and/or RL case). 
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